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The parity modulation of the ground state of a 
superconducting island is a direct consequence of 
the presence of the Cooper pair condensate pre¬ 
ferring an even number of charge carriers mm- 
The addition energy of an odd, unpaired quasi¬ 
particle equals to the superconducting gap, A, 
suppressing single electron hopping in the low 
temperature limit, ksT C A. Controlling the 
quasiparticle occupation is of fundamental impor¬ 
tance for superconducting qubits as single elec¬ 
tron tunneling results in decoherence |3l |4]. In 
particular, topological quantum computation re¬ 
lies on the parity control and readout of Majo- 
rana bound states mil- Here we present parity 
modulation for the first time of a niobium titan- 
ite nitride (NbTiN) Cooper-pair transistor cou¬ 
pled to aluminium (Al) leads. We show that this 
circuit is compatible with the magnetic field re¬ 
quirement B ^ 100 mT of inducing topological su¬ 
perconductivity in spin-orbit coupled nanowires 
[iHai. Our observed parity lifetime exceeding 1 
minute is several orders of magnitude higher than 
the required gate time of flux-controlled braid¬ 
ing of Majorana states [10 j. Our findings readily 
demonstrate that a NbTiN island can be parity- 
controlled and therefore provides a good platform 
for superconducting coherent circuits operating 
in a magnetic field. 

Experimentally, the parity modulation of a supercon¬ 
ducting island can be observed via the ground state 
charge [2], the even-odd modulation of the charge sta¬ 
bility diagram mun], or the parity dependence of the 
switching current, isw m- The interplay of the charg¬ 
ing energy Ec = 12C and the Josephson coupling 

Ej = IcTi/2e makes the Cooper-pair transistor (CPT) a 
single, gate-modulated Josephson junction USUI] with 
a 2e charge periodicity in the absence of parity switches, 
i.e. infinitely long parity lifetime, Tp. 

Recent, direct measurements of Tp mm yielded val¬ 
ues up to the millisecond regime for aluminium devices. 
Despite considerable efforts, however, no 2e periodic¬ 
ity has been reported for non-aluminium superconduc¬ 
tors USHIH], such as niobium or vanadium. Comparative 
studies of aluminium and niobium CRTs suggested that 
the elusiveness of parity effects is related to the mate- 
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FIG. 1. Device layout and parity-dependent switch¬ 
ing current, (a) Scanning electron microscope image and 
schematic cross-section of a typical Al/NbTiN/Al hybrid 
Cooper-pair transistor (CPT). The NbTiN island dimensions 
are 250 nm x 450 nm. Scale bar denotes 200 nm. The A\Ox 
oxide barriers are indicated by thick red lines, (b) Energy 
level diagram as a function of the gate charge, Ug, in the 
presence of low energy subgap states restoring 1 e periodicity. 
Gray lines denote energy levels in the absence of Josephson 
coupling, i.e. Ej = 0. Red and blue lines show energy levels 
for even and odd charge parity respectively, both for Ej — Ec. 
Parity switches occur on the timescale of Tp due to quasipar¬ 
ticle tunneling. Measured switching current histogram (c) 
and calculated Icijig) (d) in the low temperature limit. Note 
that in (d) the two possible Ic{Tig) values corresponding to 
the even and odd charge states denoted by blue and red lines 
respectively. In the measured data (c) the two branches are 
superimposed, see text. 

rial properties m, in accordance with earlier measure¬ 
ments showing that subgap quasiparticle states may ap¬ 
pear in niobium due to oxidization of the bulk material 
m- In contrast, nitridized niobium compounds, such as 
niobium titanite nitride (NbTiN), have been shown to be 
less prone to form oxides and hence are good candidates 
for parity-conserving superconducting circuits. Eurther- 
more, NbTiN forms transparent contacts with spin-orbit 
coupled semiconductor nanowires [9], and has become a 
preferred superconductor to investigate Majorana bound 
states. 

Our device features a NbTiN island sputtered onto Al 
leads (Eig. la). The tunnel barriers between the island 
and the leads are created by means of controlled in-situ 
surface oxidization of Al, resulting in amorphous AlO^, 
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FIG. 2. Characterization and temperature dependence of the parity lifetime, (a) Linecut histogram at integer Ug 
of Fig. Ic showing a bimodal distribution. We attribute the two peaks to the two parity states of the CPT (colored red and 
blue, respectively), (b) For a fast current ramp (upper panel), the histograms of the two parity states are independently 
probed showing the characteristics of data in panel (a). In the slow limit (lower panel in (b)), parity switches occur during the 
current ramp, leading to an exponential tail of the distribution (shown in green), quantifying Tp. The /bias(^) current ramp is 
represented by a dark gray line, (c) Experimental dataset in the slow limit. Note the change in the current- and timescale 
compared to panel (a). We show the exponential cutoff in green, and extract Tp = 49s (solid black line), (d) Experimental 
data at different temperatures show the temperature dependence of Tp. (e) Tp as a function of temperature for non-shielded 
device N, and shielded devices SI, S2. Eor detailed comparison, see the main text and the supplementary material. All CPTs 
exhibit an activated behavior with A 170 ... 210/iV in the high temperature limit corresponding to the gap of the aluminium 
leads. Saturation of device N without shielding and no quasiparticle traps is observed below T ^ 100 mK. Shielded devices SI 
and S2 exhibit a minigap-activated behavior A* 20 jaV in the low temperature limit. The fitted Tp(T) function is defined in 
the main text. Star symbol shows Tp = 15s at T = 12mK extracted from parity distilled data for device S2 (see Eig. 3d). 


barriers m- 

We extract a charging energy Ec ~ 50/ieV from the 
measured charge stability diagram. For different devices, 
we estimate Ej « 30... 50 /ieV from the superconduct¬ 
ing gaps and normal state resistances of the junctions m 
assuming equal resistances for the two tunnel barriers. A 
detailed list of parameters and characterization methods 
are presented in the supplementary material. Our devices 
are in the intermediate coupling regime with Ej r\j Ee, 
where the energy diagram (Fig. lb) and the critical cur¬ 
rent (Fig. Id) are sensitive to the charge parity. It is im¬ 
portant to note that our CPTs are in the optimal regime 
to establish flux-controlled braiding of Majorana bound 
states with Ej ^ Eq k^E uni and hence a useful plat¬ 
form to establish the parity lifetime for Majorana circuits 

m 

We model the CPT as a two level system which can 
exist in either parity state (red and blue bands in Fig. 1, 
respectively), and switches state on the timescale of Tp 
|23| due to quasiparticle tunneling. We collect the switch¬ 
ing current histograms by repetitively sweeping the bias 


current from zero (non-dissipative state) to beyond the 
switching current. Here, in the resistive state, quasiparti¬ 
cle tunneling causes a random reinitialization of the par¬ 
ity state of the CPT for the next measurement. This re¬ 
sults in the apparent le periodicity in Fig. Ic. Neverthe¬ 
less, as long as the parity remains constant during each 
sweep, we expect to find the two branches as a bimodal 
histogram, as we indeed observe in Fig. 2a. In these mea¬ 
surements, the current ramp time is much shorter than 
the parity lifetime, Tp {fast measurement limit). 

We quantify Tp in the slow measurement limit. In 
this regime parity switches occur during the current 
ramp (Fig. 2b lower panels) such that reaching the up¬ 
per branch (depicted as blue in Fig. 2a and 2b) be¬ 
comes exponentially suppressed (Fig. 2c). The expo¬ 
nential tail represents parity switches during the cur¬ 
rent bias ramp, resulting in an observable decay of the 
upper branch (depicted as green in Fig. 2b and 2c), 
Pu{t) = Pw(0) exp(—t/Tp). Thus, from the decay of the 
histogram (black solid line in Fig. 2c), we can directly 
obtain Tp. 
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The observed is a result of single electron tunneling 
events through the junctions of the CPT. In the zero 
voltage state, we can link Tp to the subgap density of 
states on the island and the quasiparticle density in the 
leads [24] (for details see the supplementary material): 

^-1 _ 2nqp Psubgap 
^ e‘^RNpA\ PNbTiN ’ 

where Rn is the normal state resistance of the CPT. It is 
instructive to note that the parity lifetime is determined 
by the quasiparticle density in the leads [22] and the phe¬ 
nomenological Dynes parameter [25] of the island mate¬ 
rial. Assuming a thermal nqp(T) oc VT exp(— 
in the leads, we find Aai = 170... 210 peY (Fig. 2e) 
for temperatures exceeding 120 mK, in good agreement 
with the superconducting gap of the aluminium leads 
extracted from the charge stability diagram. We can 
therefore attribute the observed parity lifetime for T > 
120 mK to the thermal quasiparticle population in the 
leads. 

For device N, however, we find a saturated Tp = 9.5 ms 
in the low temperature limit, a common observation in 
superconducting qubits [4] and hybrid single electron 
transistors [24] signifying the presence of non-thermal 
quasiparticle excitations. 

We further improve the low temperature parity lifetime 
by introducing microwave-tight shielding coated with in¬ 
frared absorber painting and Ti/Au quasiparticle traps 
for devices SI and S2. These additions result in a non- 
saturated behaviour of r^, and we observe a minigap ac¬ 
tivated behaviour with A* ^ 20 peY for both devices. It 
is to be stressed that this observation signals that the 
effective quasiparticle temperature of the CPT follows 
the bath temperature down to the 10 mK regime. We 
find Tp = 49 s at T = 12 mK for device S2. To put 
this number into context, we note that the Josephson 
frequency fj = Ej/h^ 10 GHz and thus a single quasi¬ 
particle event occurs only once for every Tpfj ^ 10^^ 
Cooper pairs tunneling through the junctions. This sig¬ 
nifies the low probability of parity switches for an open 
device with Ej « Ec required for flux-tunable Majorana 
braiding schemes. 

Thus far, we started each switching current measure¬ 
ment from an unknown parity state because of the ran¬ 
dom reinitialization in the dissipative state of the CPT 
for /bias > 7sw In order to reproducibly select the 
same parity state, we employ a parity distillation pro¬ 
tocol (Fig. 3) where, by selecting a single parity state 
without switching to the resistive state, we ensure that 
the parity remains well defined for the subsequent mea¬ 
surement. This protocol indeed results in a 2e periodic 
switching current pattern (Fig. 3c) which is observed for 
the first time for a non-aluminium CPT. 

We quantify the effectiveness of the parity distillation 



FIG. 3. Parity state distillation, (a) Schematic cur¬ 
rent bias and gate voltage waveforms applied to the device. 
First, current bias pulse with an amplitude between the two 
branches is applied (1). At the same time the measured volt¬ 
age Vsd{t) is recorded to initialize the parity state. Then the 
gate voltage is ramped (2) and finally, the switching current is 
recorded (3). (b) The schematic representation of the applied 
waveforms overlayed on the parity-dependent critical current, 
(c) Parity distilled switching current histogram exhibiting 2e 
periodicity with twait = 100 ms. (d) Parity distillation fidelity 
as a function of twait, see main text. Solid red line denotes the 
fit ~ exp(—twait/Tp) with Tp = 15 s. All data was recorded on 
device S2. 


by defining the fidelity as: 


Tp = 


Pu,f Pu,i 
1 - Pu,i 


( 2 ) 


where Puj is the conditional probability of the upper 
branch in the final step (3), and Pu,i is the initial prob¬ 
ability. We note that the above expression is valid for 
an arbitrary 0 < Pu,i < 1 value set by the average quasi¬ 
particle occupation of the CPT. For device S2 we find 
Ep = 0.88 ± 0.05 for t^ait = 100 ms demonstrating the 
high degree of parity distillation. By changing twait be¬ 
tween the parity initialization (1) and measurement (3), 
we acquire an independent measurement Tp = 15 s for 
device S2 at T = 12mK (Fig. 3d). 

Finally, we investigate the evolution of Tp{B) in differ¬ 
ent magnetic field directions. In parallel magnetic field, 
we observe a gradual decrease of Tp. The onset of the 
steep decay at H = 110 mT (yellow arrow in Fig. 4a) is 
in agreement with the condition for vortex penetration 
through a mesoscopic superconducting island [26] with 
<1> ^ 1.14>o ^ 4>o. This result underlines the significance 
of the sample geometry for magnetic field-enabled CPTs 
to avoid vortex formation. However, our device exhibits 
Tp > 10 ms in H|| > 100 mT, required to induce Majorana 
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FIG. 4. Influence of the magnetic fleld on the CPT. 

(a) The parity lifetime as a function of the in-plane field, 
pointing as shown in the inset. The shaded region depicts 
the condition for inducing Majorana bound states (MBS) in 
spin-orbit coupled semiconductor wires [3E1- The change in 
slope dil B ^ 110 mT is denoted by the yellow arrow and the 
dash line serves as a guide to the eye. (b) The parity lifetime 
as a function of perpendicular field, B^. The dash line shows 
the overall trend. The inset shows the electron microscope 
image of the device with the white scale bar denoting 5 /xm. 


bound states [7H9]. 

In contrast, applying a perpendicular field first results 
in an increase of reaching a maximum value Tp = 94 s 
at B^ ^10 mT, before dropping at higher magnetic fields 
(Fig. 4b). Making use of the relation between the lower 
critical field and the stripe width, Bci{w) ^ [27] 

we attribute the initial increase to vortex formation and 
hence more effective quasiparticle trapping in the wide 
lead sections {w ^ 2 /xm) of the device (see wide light 
gray sections in Fig. 4b). Upon reaching B ^ 10 mT, 
the vortex phase becomes stable in the close vicinity of 
the CPT {w ~ 250 nm) causing a gradual decrease of Tp. 
We note that due to ^ai ~ 100 nm ^ w for our devices, 
we cannot make a quantitative comparison between our 


measurement data and theoretical expressions of criti¬ 
cal magnetic field of thin stripes. We however conclude 
that vortices induced by a perpendicular magnetic field 
can increase the efficiency of quasiparticle traps, but the 
formation of a vortex phase in the near vicinity of the 
CPT enhances quasiparticle transport in agreement with 
earlier observations [28l [29] . 

In conclusion, we fabricated and characterized NbTiN 
Cooper-pair transistors showing parity effects for a non¬ 
aluminium superconductor for the first time. We char¬ 
acterize the parity lifetime by evaluating the switching 
current histograms and find values exceeding 1 minute at 
T = 12 mK. Our devices are in the regime of Ej ^ Ec^ 
ksT ideal for Majorana braiding schemes. Furthermore, 
we demonstrate charge parity distillation to reproducibly 
initialize the island in a given parity state. We find par¬ 
ity lifetimes in excess of 10 ms in external magnetic fields 
up to 150 mT and showed the importance of sample ge¬ 
ometry for magnetic field-enabled operation required for 
inducing Majorana bound states. 
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Supplementary online material 
1 minute parity lifetime of a NbTiN Cooper-pair box 

DEVICE FABRICATION 

The Cooper-pair transistors (CPTs) were fabricated using electron beam lithography and thin film deposition 
starting with a p++ doped silicon wafer with a 285 nm thick thermally grown Si02 surface layer. First, aluminium 
leads were defined and evaporated in a high-vacuum chamber (pbase ^ 10“^ Torr) at a rate of 0.2 nm/s with a thickness 
of 30... 35 nm. Subsequently, the mask for the NbTiN island was defined in a second lithography step. The sample 
was loaded into an ultra high vacuum AJA International ATC 1800 sputtering system (pbase ^ 10“^ Torr), where first 
a ^ 5 nm A1 layer was removed by means of argon plasma etching at p = 3 mTorr. This step was followed by in-situ 
oxidization to create AlO^, tunnel barriers. Without breaking vacuum, the NbTiN island was then sputtered with a 
layer thickness of 70 ... 100 nm. 

We used a Nbo.yTio.s target with a diameter of 3". Reactive sputtering resulting in nitridized NbTiN films was 
performed in an Ar/N process gas with a typical 10at% nitrogen content at a pressure of 3 mTorr using a DC 
magnetron source. During deposition, a —45 V bias was applied to the sample with respect to the target. The critical 
temperature of the superconducting transition temperature of thin films with a layer thickness of 100 nm was measured 
to be 14.1 K in zero magnetic field. 

For the shielded samples SI and S2, quasiparticle traps were fabricated by first cleaning the A1 surface by means 
of argon plasma milling at p = 0.2 mTorr and in-situ evaporation of 25 nm Ti and 100 nm An films at a base pressure 
of ^ 10“^ Torr (see Fig. Sic and d). 

Care was taken to remove resist mask residues after each electron beam writing step using a remote oxygen plasma 
etch with a pressure of 1 mbar. 


device 

island size 
(nm X nm) 

junction size 
(nm X nm) 

NbTiN thickness 
(nm) 

traps 

oxygen exposure 
(Torrxs) 

Rn 

(kn) 

N 

500 X 500 

200 X 200 

70 

no 

7400 

58 

SI 

450 X 200 

200 X 200 

100 

yes 

150 

125 

S2 

450 X 250 

200 X 250 

100 

yes 

150 

66 


TABLE SI. Device fabrication parameters of the CPTs discussed in the main text. 


MEASUREMENT SETUP 

The measurements were performed in a Leiden Cryogenics CF-1200 dry dilution refrigerator with a base temperature 
of 12 mK equipped with Cu/Ni shielded twisted pair cables thermally anchored at all stages of the refrigerator. 

Fig. S2a shows the schematics of the measurement. The current bias and gate voltage were applied through battery 
operated and optically isolated sources in order to reduce line interference. Similarly, the first stage of the Vgd 
amplifier was isolated from the commercial readout electronics. Filtering of the measurement lines was achieved by 
room temperature LC tt filters with a cutoff frequency of ^ 100 MHz followed by a sequence of a high frequency 
copper powder filter m and a two-pole RC filter with a nominal /-mb = 50 kHz, both thermally anchored to the 
12 mK stage. 

Special care was taken to avoid stray microwave radiation by using an outer and an inner copper shield enclosing the 
device. The inner surface of both shields was treated with commercially available Aeroglaze Z306 paint m absorbing 
far infrared stray radiation [32]. We note that the inner shield was not present for device N. 

A typical DC V — I trace of device S2 at 12 mK temperature is presented in Fig. S2b exhibiting a sharp transition 
between the dissipationless and the resistive state at the switching current, /gw We observe a retrapping current 
Ir <C Isw characteristic to unshunted Josephson junctions in the low temperature limit [33|. The additional features 
in the resistive state are consistent with Fiske steps [34] . 

The switching current histograms were acquired using a Rigol DG4062 arbitrary waveform generator controlling the 
isolated current bias source with a triangle wave signal resulting in a dl/dt current ramp. A finite voltage response 
above the preset V^ei ^ 10/iV triggers the recording of Aw (Fig- S2c). We note that the delay of the low pass filters 
were accounted for on the basis of circuit simulations. Subsequent Aw measurements were taken with setting zero Aias 
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FIG. SI. Scanning electron microscope images of the CPTs. (a) Non-shielded device N without quasiparticle traps, 
(b), (c), (d) Shielded device S2 featuring quasiparticle traps. Colour legend: dark red: Si/Si02 substrate; light blue: NbTiN 
island, green: aluminium leads and yellow: normal metal (Ti/Au) gate and quasiparticle traps. Scale bars denote 500nm (a), 
1000 nm (b), 5/xm (c) and 200 nm (d), respectively. 


for approximately 100 ms in between to avoid overheating effects. We did not observe a difference in the switching 
current histograms taken with waiting times in the range of 20 ms and 10 s. 


BASIC CHARACTERIZATION OF THE CPT 

In Fig. S3a we present a typical charge stability diagram of the CPT. First, we establish the superconducting gap 
of the aluminium lead and NbTiN island by finding the onset of quasiparticle transport (Fig. S3b and c). 


device 

^Al 

(MeV) 

^NbTiN 

Rn 

(kn) 

Ej 

Ec 

(MeV) 

Ej/Ec 

calc. 

Ej/Ec 

^SW 

N 

210 

1390 

58 

54 

50 

1.08 

1.25 

SI 

218 

1420 

125 

21 

62 

0.34 

0.32 

S2 

220 

1300 

66 

48 

49 

0.98 

1.16 


TABLE S2. Device transport parameters. 

It is important to note that we find a finite subgap conductance at eVsd < 2ANbTiN which is consistent with 
the presence of the subgap quasiparticle states justifying the analysis leading to equation (1) in the main text. 
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FIG. S2. Measurement electronic setup and typical waveforms, (a) Schematics of the measurement, (b) Typical 
experimental V — I trace exhibiting a well-defined switching current, /sw (c) Current ramp (red) with a slope of dl/dt probing 
Isw which is recorded when the measured Vsd (blue) reaches Vref ~ 10 /iV. 



FIG. S3. Basic characterization of the CPT. (a) Stability diagram: log |/| as a function of Vsd and Vg. We estimate Ec 
based on the characteristic resonances (see text). (b),(c) I — V (black) and dljdV (red) line traces on the scale of Aai and 
ANbTiN, respectively. All data was taken for device SI. 


Furthermore, we observe reduced ANbTiN ~ 1.3... 1.4meV values compared to that of bulk films (A > 2meV) [35] . 
which we attribute to the chemical interaction between the AlO^, tunnel barrier and the NbTiN film. Indeed, it was 
shown earlier that the critical temperature of Nb films is particularly very sensitive to contamination with oxygen 
|36l|37]. However, the nitridized NbTiN compound is presumably less prone to oxidization [38] . 

We evaluate the Josephson coupling for a single tunnel junction [39] : 


Ej 


T__^alk ( li ( 

2e^RN/2 ly V ^NbTiN / I 


(SI) 


with K{x) being the complete elliptic integral of the first kind. This expression is valid in the zero temperature limit. 
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assuming symmetric tunnel junctions of the resistance of Rn/2. 

We estimate the charging energy, En = e^/2(7, based on the periodicity of characteristic resonances visible for 
eysd<2AAi (Fig. S3a) [40l|4l]. 

Alternatively, we can estimate the EjjE^ ratio based on the modulation of /gw as a function of the gate charge, rig 
(last column of Table S2) [42]. We find a good agreement between the Ej jE^ values acquired by the two independent 
methods. 


EVALUATION OF THE PARITY LIFETIME 



FIG. S4. Measured switching current histograms at different dl/dt ramp rates. Colour legend: red: dl/dt = lOOnA/s, 
black: dl/dt = 400nA/s. The data was acquired on a CPT similar to device N. 


We checked the robustness of the extracted parity lifetime against changing the current ramp rate. Typical datasets 
are shown in Fig. S4, giving = 1.2 ms and Tp = 0.98 ms for dl/dt = lOOnA/s (red) and dl/dt = 400nA/s 
(black), respectively. We estimate the typical uncertainty to be 15%, concluding that Tp does not depend on dl/dt 
which validates the analysis in the main text. However, we do not discuss here the intrinsic peak shapes of the 
bimodal switching current histogram. Since fast gate charge noise influences the measured distribution [43], we 
cannot distinguish between thermally activated m and macroscopic quantum tunneling behaviour [331 [45]. 

We now turn to the temperature dependence of Tp. In order to get equation (1) in the main text, we assume the 
following: 

i) The superconducting gap of the island (ANbTiN ^ 1.3 meV) is much higher than the effective thermal energy 
describing the quasiparticle population and the gap of the leads (Aai ~ 200 /reV). 

ii) The density of states in the leads is BCS-type: piead(F^) = PA\ ^\E\/ for \E\ > Aai and zero otherwise. 

iii) There is a constant, finite subgap density of states Psubgap for energies below ANbTiN in the island. 

iv) the energy dependence of the single electron tunnel probability is negligible over the energy range of ^ ANbTiN, 
meaning that the tunnel barrier is much higher than ANbTiN- 

v) The tunnel barriers are identical, each characterized by half the normal state resistance of the full device, R^. 

Considering only single electron tunneling and zero voltage bias across the tunnel barriers, following m, we get 
the quasiparticle tunnel rate: 


= 


Psubgap 


e‘^{RN/2)pA\ PNbTiN Jo 


f 


Pxe.o.{E)f(T,E)dE = 


277/qp Psubgap 
C^RnPAX PNbTiN- 


(S2) 
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Assigning an effective temperature to the quasiparticle population in the leads, we find: 

. _ aai 

nqp{T)=2 / p\Qg,(^{E)f{T,E)dE = pA\\/27ikBTAA\e (S3) 

Jo 

in the limit of ksT Aai- For temperatures exceeding 100 mK, we assume that the quasiparticle population is in 
thermal equilibrium, and therefore the lattice temperature is equivalent to the effective quasiparticle temperature: 
Tp^{T) oc y/r exp(-AAi/fcBT). 

We verify this picture by fitting the observed parity lifetimes with Aai as a free parameter, and find values ranging 
170... 210 /ieV for different devices in good correspondence with the gap determined by voltage bias spectroscopy 
(Aai in Table S2). 

Notably, the ratio Psubgap/PNbTiN is the Dynes parameter [47] of the island material, characterized to be ^ 10“^ 
based on measurements of highly resistive single junctions. With this value and using pA\ = 1-45 x 10^^ |46| . 

we get riqp ~ 3 for device N based on the observed parity lifetime of 9.5 ms in the low temperature limit. 

We now comment on the observed A* 20 /reV activation energy observed for devices SI and S2. We estimate the 
maximal even-odd energy difference to be 6E ^ 20... 30 /reV based on Ej ^ Ec ^ 50 /reV [48] which is in range of 
the experimentally observed A*. Similar, activated behaviour of the parity lifetime scaling as ^ ex.p{SE/kBT) was 
reported earlier |49| . 

Providing another possible explanation, we note that a grain size of ~ 50 nm can lead to a level spacing of the 
order of 10 /ieV which can influence single electron transport and hence Tp if the grains are weakly coupled, i. e. for 
disordered superconducting films m- Disorder-induced fluctuations may also explain the broadening of the coherence 
peaks (Fig. S3c) [571 [5^ . 


SUPERCONDUCTING THIN FILM CHARACTERIZATION AND MAGNETIC FIELD DEPENDENCE 


Next, we consider the properties of superconducting stripes with layer thickness d, and a width w to find the London 
penetration depth Xl and the coherence length We characterize the upper critical field in the parallel {Bc\\) and 
perpendicular geometry {Bc±) based on the dl/dV traces of the tunnel junctions of the CPT. In addition, we measure 
the normal state resistivity of the films that gives an estimate for the mean free path, Im m- 

First, we establish the length scales of the island material, NbTiN. We find films superconducting at B± = 9T 
which leads to an upper limit of ^NbTiN < b nm following m-- 


Bc± = 


^0 

27r^2' 


(S4) 


The penetration depth can be estimated using the normal state resistivity of p = 98 pQcm and the critical temperature 
of T = 14.1 K using the following semi-empirical formula [35] : 


ANbTiN = Y [K] ^ ^ ^ 

Next we estimate length scales of the A1 leads based on the electronic transport through the CPT. Typical thin A1 
films are type-II superconductors in the dirty limit {Im < Co) with a reduced coherence length of C ~ 0.85 a/Co^ and 
with a London penetration depth of A ^ Ao^/Co/^m^ where Co ~ 1500 nm and Aq ~ 16 nm are the bulk values [5l]. For 
our films, we estimate Im ~ 8.5 nm based on the resistivity of p = 4.3 pUcm m- From the stability diagram of the 
devices, we extract upper critical fields of Bc 2 ,i. = 36.4 ± 4 mT (Fig. S5d) and 5c2,|| = 320±10mT (Fig. S5b) leading 
to a coherence length of Cai = 96 nm and Aai = 230 nm which enables vortex formation in the aluminium leads in 
perpendicular magnetic field. 

It is important to observe that Cnotin ^ d ~ 100 nm enables vortex formation for an in-plane geometry in the 
NbTiN island. We find a characteristic suppression of Tp at B\\ = 70 mT for device S2 (orange dots in Fig. S5a) and 
at B\\ = 110mT for device SI (cyan dots in Fig. S5a). Considering the effective cross-sectional areas (see Table SI for 
dimensions), we find ^ 1.5<Fo and 4> « 1.14>o for S2 and SI, respectively, which is in qualitative agreement with the 
threshold of a single vortex formation in a mesoscopic island [55H57]. We also note that 5c2,|| of iFe leads (Fig. S5b) 
does not depend on the direction of B \\, therefore the different evolution of Tp can only be explained by the different 
alignment of B\\ with respect to the NbTiN islands. 
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FIG. S5. Parity lifetime and superconducting gap in magnetic field, (a) Measured parity lifetime as a function of 
in-plane magnetic field. The inset shows the relative orientation of the field for the two devices, (b) The superconducting gap 
of the A1 leads as a function of in-plane magnetic field. (c),(d) The parity lifetime and the superconducting gap as a function 
of the perpendicular magnetic field. Colour legend: cyan: device SI; orange: device S2. The scale bar in the inset of panel (c) 
denotes 5 /xm. 


In a perpendicular geometry, the vortex phase is stable in a thin stripe above the magnetic field 

Bci,±{w)=a^, (S6) 

where a is a model-dependent prefactor [58H6Q] of the order of unity. We reproducibly find the same non-monotonic 
behaviour of for devices SI and S2 with the maximum at ^ 10... 13 mT, which is in range of for 

w = 200 ... 250 nm, the width of the A1 leads near the island. 
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